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mis 
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影响 水 文 和 气候 变化 的 基本 参数 ,研究 土壤 含水 量 分 布 ， 对 气候 变化 、 水 资源 分 布 、 农 


是 
作物 估 产 等 有 着 重要 的 现实 意义 和 科学 价值 。 本 文 以 2015 年 6 月 21 日 的 Sentinel-1A( 哨 兵 1 号 ) 双 极 化 合成 
孔径 雷达 影像 为 基础 ， 结 合同 时 段 辅助 光学 影像 Landsat 8， 对 北 安 一 黑河 高 速 沿 线 地 区 不 同 植被 覆盖 程度 下 
复杂 地 表土 壤 含 水 量 进行 反 演 研究 ， 探讨 不 同 极 化 组 合 方式 在 不 同 土地 利用 方式 下 的 土壤 水 分 含量 反 演 结 


果 。 结 果 表 明 : VH 极 化 及 


VH 5 fi Ej 7x € NDVI(Normalized Difference Vegetation Index) 组 合 反 演 精度 分 别 为 


52.1% 和 53.6%， 整 体 效 果 并 不 理想 。VV 极 化 (VV Polarizatiom) 图 像 和 双 极 化 VV/VH(VH Polarization) 组 合 在 
裸露 和 低 植被 地 区 反 演 更 具有 优势 ， 其 精度 分 别 为 75.4% 和 59.5%， 而 在 高 植被 覆盖 度 地 区 并 不 适用 。VH 极 


化 反 演 结果 中 耕地 土壤 含水 量 比 实 际 值 低 9.37%, VV 极 化 在 低 植被 区 域 土壤 含水 量 比 实际 值 低 10.45%, 在 灌 


木 及 耕地 地 区 VV/VH 反 演 结 


果 精 度 比 单 极 化 及 其 组 合 反 演 结果 低 ， 最 高 精度 模型 的 反 演 是 VV 结合 NDVI. 


VV 与 辅助 变量 NDVI 结合 能 综合 反映 复杂 地 表 环 境 下 土壤 含水 量 ， 其 精度 达 84%， 标 准 误差 RMSE 为 2.07, 
比 VV 极 化 反 演 精度 提高 8.8%, RMSE H VV 极 化 降低 2.704. VV 与 辅助 变量 NDVI 组 合 方式 在 中 等 植被 履 
盖 地 区 土壤 含水 量 反 演 更 具有 优势 ， 并 能 够 更 好 地 发 挥 哨兵 1 号 C 波段 合成 孔径 雷达 在 土壤 水 分 研究 中 的 法 


力 与 有 效 性 。 
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Abstract: Soil water content is one of the basic parameters that affect hydrological variability and climate change. It has 


important practical significance and scientific value for climate change, water resources and estimation of crop yield to study 


the distribution of soil water content. To probe new ways of soil water content retrieval in complex vegetation coverage area, 
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this study analyzed soil water contents of complex surfaces with various degrees of vegetation cover along Bei'an-Heihe 
Expressway using images from Sentinel-1A dual-polarization Synthetic Aperture Radar (SAR) for 21 June 2015. Also Landsat 
8 images were integrated as assisted optical image for the same satellite transit time. Then the results of the inversion of soil 
water content under different land use types and polarization combinations were discussed. Backscattering coefficients of 
different polarization modes were extracted using the water cloud model. Support vector regression algorithm was used to 
estimate surface soil water content based on the soil inversion parameters. The applicability of different polarizations in the 
retrieval of soil water content on complex surface was also discussed. The results showed that VH polarization retrieval 
accuracy was 52.1196, while combined VH polarization with normalized difference vegetation index (NDIV) retrieval accuracy 
was only 53.696. This was not satisfactory for the vegetation zone. VV polarization and dual polarization ratio of VV/VH 
images were very sensitive to bare land and low vegetation cover land, for which retrieval accuracies were respectively 75.4% 
and 59.5%. These methods were, however, not applicable in areas with moderate or high vegetation cover. The results of VH 
polarization inversion for arable lands soil water content was 9.37% lower than the measured value. Also the inversion value of 
VV polarization for areas with low bush was 10.45% lower than the measured value. The inversion results for dual polarization 
ratio of VV/VH in shrub and arable lands were not as good as the inversion results for single polarization. For the various 
combinations, the inversion with the highest precision model was that for the combination of VV with NDVI. In summary, the 
combination of VV and auxiliary variable NDVI comprehensively reflected soil water content in complex surface 
environments. The goodness of fit (R°) of VV polarization combined with NDVI was 84% and the calculated root mean 
squared error was 2.07. In comparison with VV polarization, the retrieval accuracy improved by 8.8% and the calculated root 
mean square error decreased by 2.704. The combination of VV polarization with NDIV had more advantages for the inversion 
of soil water content for the regions with middle vegetation cover. The application of combined VV polarization with NDIV 
increased the potential and effectiveness of Sentinel-1A c-band synthetic aperture radar in areal study of soil water content. 


Keywords: Sentinel-1A; Soil water content; Backscattering coefficient; Polarization; Support vector regression 


不 易于 使 用 。 研 究 表明 当 微波 遥感 数据 协同 光学 植 
被 指数 的 方式 "1 反 演 土壤 合 水 量 时 ,用 光学 数据 计 
算 植被 生物 物理 参数 , 能 更 好 地 补充 植被 的 影响 。 
国内 外 研究 虽 在 不 同 植被 覆盖 区 域 土壤 合 水 量 反 演 


土壤 合 水 量 (soil water content，SWC) 是 影响 水 
文 和 和 气候 变化 的 基本 参数 ， 它 不 仅 是 联系 地 表 水 与 地 
下 水 的 纽带 ,也 是 研究 地 表 能 量 交换 的 基本 要 素 M, 在 
气候 学 中 , 土壤 含水 量 的 变化 能 影响 其 本 身 的 水 热 过 


程 ,使 地 表 参 数 发 生变 化 从 而 影响 到 气候 的 变化 中 在 
生态 学 中 ,土壤 含水 量 在 土壤 -植被 -大 气 连续 体 物 
质 与 能 量 转化 中 也 发 挥 着 重要 作用 FJ， 直接 影响 生 
态 环境 的 组 成 。 因 此 ， 研 究 土壤 含水 量 分 布 ， 对 气候 
变化 、 水 资源 分 布 、 生 态 退 化 过 程 以 及 农作物 估 产 
等 有 着 重要 的 现实 意义 和 科学 价值 。 

目前 用 于 土壤 含水 量 反 演 的 遥感 主要 有 3 类 : 
可 见 光 外 、 热 红外 双 感 中 和 微波 么 感人 中 。 其 中 ,微波 
遥感 具有 对 土壤 水 分 的 高 敏感 性 ， 以 及 测量 不 受 云 
层 和 黑夜 影响 的 能 力 ， 在 土壤 含水 量 估 测 中 应 用 非 
常 广泛 5。 近 年来， 在 裸露 地 表 通 常 使 用 一 些 物 理 经 
验 模 型 (如 Oh RÆ, Dubois 模型 中 进行 土壤 合 水 
量 估 测 ， 但 需要 现场 测量 进行 校准 ; 在 低 植 被 覆盖 
度 地 区 ， 通 常 选 择 将 半 经 验 模型 (如 Shi 模型 路、 水 
云 模 型 09 等) 应 用 到 水 分 反 演 中 ,植被 区 土壤 含水 量 
仍 被 过 低估 计 ; 对 于 植被 高 覆盖 度 地 区 ， 植 被 冠 层 
双向 散射 衰减 了 土壤 后 向 散射 的 影响 ， 需 要 采用 密 
歇 根 微波 植被 散射 模型 (Michigan Microwave Canopy 
Scattering, MIMICS)05 及 其 改进 模型 (如 Bi-MIMICS 
模型 中 计算 裸 土 后 向 散射 系数 ， 由 于 参数 太 过 繁琐 


中 有 不 同 的 方法 , 但 总 体 缺 乏 针 对 复杂 地 表 较 为 综 
合 的 简单 方式 。 

本 研究 以 北 黑 高 速 路 段 治 线 区 域 为 案例 区 ， 由 
于 受 全 球 气候 变化 的 影响 ,这 一 地 区 属于 古 冰川 沉 
积 残留 物 的 岛 状 冻 土 正 发 生 退 化 ， 地 表 环 境 破碎 且 
复杂 。 采 用 水 云 模型 对 sentinel-1A 合成 孔径 雷达 
(synthetic aperture radar，SAR) 影 像 提 取 其 不 同 极 化 
条 件 下 的 土壤 后 向 散射 系数 ， 利 用 光学 影像 作为 畏 
助 数据 提取 复杂 地 表 的 植被 归 一 化 指数 ， 对 上 比 支 持 
向 量 回 归 算 法 下 不 同 极 化 与 光学 参数 之 间 的 反 演 组 
合 精 度 并 分 析 适 用 性 情况 ,为 不 同 植被 覆盖 地 区 土 
壤 合 水 量 反 演 研究 提供 新 的 思路 ， 同 时 为 在 高 纬度 
多 年 冻 土地 区 开展 土壤 合 水 量 监测 提供 科学 依据 。 
1 数据 来 源 与 研究 方法 
14 ”研究 区 概况 

研究 区 为 北 安 至 黑河 ( 北 黑 ) 高 速 公 路 治 线 
(127?17'31"E-127?31'24"E, 49?30'57"N-49?41'50"N), 
位 于 黑龙 江 省 小 兴安 崔 的 西北 部 (图 1)。 处 于 高 纬度 
多 年 冻 土 区 与 季节 性 冻 土 过 渡 地 带 ， 具 有 典型 冰 缘 
地 貌 ， 岛 状 冻 土 活动 层 厚 度 约 为 1.8 m。 山 谷 、 低 洼 
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1 研究 区 地 理 位 置 (a) 及 样本 点 分 布 (b) 


Fig. 1 


河流 阶地 、 阴 坡 和 湿地 为 岛 状 多 年 冻 土 的 发 展 提供 了 
合适 的 条 件 ， 岛 状 冻 土 在 S 月 深度 最 大 , 达 2.26~2.67 m, 
融化 期 为 6 月 至 9 月 ,冻结 期 为 10 月 至 翌年 $ 月 , È 
要 分 布 集中 在 山谷 和 阴 坡 地 带 。 地 形 属于 河谷 地 段 ， 
主要 分 布 有 丘陵 、 漫 岗 、 耕 地 、 沼 泽 地 等 类 型 。 气 候 
属 山区 高 寒 地 带 ， 年 平均 气温 -2~1 C, 年 降雨 量 
500—600 mm, 年 日 照 2 500 h, 无 霜 期 约 90 do 红色 
边疆 农场 二 分 场 在 该 研究 区 内 ， 耕 地 分 布 在 坡地 及 
河谷 两 岸 ， 农 作物 以 小 麦 (Tyiticum aestivum) H Èo 


Location of the study area (a) and the distribution of sample points (b) 


主要 植被 类 型 是 森林 植被 ( 针 阔 温 交 林 ) 和 草 旬 植被 
ELFA, ERES) ARA A A tE (Betula 
platyphylla), 3&b(Tilia tuan), YEBW(Xylosma racemosum), 
落叶 松 (Larix gmelinii)、 杨 树 (Populus sp.)o 
12 ”数据 来 源 及 处 理 

土壤 含水 量 反 演 过 程 涉及 的 主要 数据 包括 SAR 
影像 、 MODIS 数据 产品 、Landsat 8 OLI 数据 、DEM 
数据 及 地 表 覆 盖 类 型 图 ,数据 分 辩 率 、 数 据 来 源 及 
数据 的 应 用 见 表 1。 


R1 北 黑 高 速 路 段 地 区 土壤 含水 量 反 演 数 据 来 源 


Table 1 Dataset list of soil water content inversion along the Bei'an-Heihe Expressway 
数据 E 用 途 来 源 
Data Resolution (m) Application Data source 
Sentinel-1A SAR 影像 5x90 提取 后 向 散射 系数 欧洲 空间 局 
Sentinel-1A SAR image Extracting the backscattering coefficients European Space Agency 
MODI IAT 地 表 温 度 产品 oTo a esci NR 
; 1 000 Calculating the empirical values of sample United States National Aeronautics and 
MODIIAI daily LST product : mn ; 
points Space Administration 
数字 高 程 模型 d SAR 地 理 编码 国际 科学 数据 服务 平台 
DEM Geocode of SAR images International Scientific Data Service 
Landsat 8 OLI 数据 30 计算 归 一 化 植被 指数 美国 地 质 勘 探 局 
Landsat 8 OLI image Calculating the NDVI United States Geological Survey 
地 表 覆 盖 类 型 图 30 登 加 分 析 全 球 地 表 覆 盖 下 载 中 心 


Land cover map 


1.2.1 Sentinel-1A 影像 及 处 理 方 法 

欧 空 局 最 新 发 射 的 Sentinel-1A 卫星 ,采用 多 种 
成 像 模式 在 土壤 含水 量 动态 监测 领域 发 挥 不 可 忽视 
的 作用 。 本 研究 所 用 数据 是 干涉 宽 模式 (Interferometric 
Wide, IW) 下 的 S1 TOPS-mode SLC 数据 ， 能 解决 宽 
晶 合 成 孔径 雷达 成 像 出 现 的 scalloping 效应 并 增强 
成 像 辐射 性 能 ,在 提取 后 向 散射 系数 中 能 够 更 精 
确 。 获 取 时 间 为 2015 年 6 H 21 A, 波段 为 C 波段 


Overlay analysis 


Global Land Cover Download Center 


(0.055 466 m),， 幅 宽 为 250 km, 分 辩 率 为 5 mx20 m, 
采用 双 极 化 方式 (VH 极 化 、VV 极 化 )。Sentinel-1A 
在 研究 区 上 空 为 右 侧 视 的 降 轨 成 像 ， 影 像 为 左右 倒 
E. 采用 GAMMA 软件 对 数据 进行 预 处 理 和 地 理 编 
B, 以 及 提取 后 向 散射 系数 ， 过 程 如 下 : 

1) 预 处 理 。 干涉 SAR 处 理 器 (ISP) 模 块 对 数据 进 
行 预 处 理 。 依 据 原始 图 像 中 的 数据 文件 (TIFF file). 
元 数据 文件 XML file)、 定 标 文件 (XML file), Ez 
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文件 (XML file) 对 相位 偏 移 和 干涉 参数 偏差 进行 辐 
射 定 标 及 去 噪 处 理 (图 2a) 其 次 是 多 视图 像 拼 接 ， 需 
要 设置 图 像 的 方向 项 及 距离 项 的 视 数 ， 将 单 视 复数 
图 生成 多 视图 像 (图 2b)。 

2) 地 理 编码 。 运 用 差分 干涉 和 地 理 编码 (DIFF & 
GEO) 模 块 对 多 视图 像 进行 地 理 编码 ， 入 射 角 为 
33.548°%。 包 括 生 成 初始 查找 表 、 影 像 配 准 、 精 化 查 
找 表 等 过 程 ， 最 后 获得 地 理 坐 标 下 的 后 向 散射 图 像 
及 SAR 坐标 下 的 DEM。 利 用 SAR 脉冲 回 波 时 延 余 
距 和 多 普 勒 中 心 频率 ， 通 过 和 斜 距 多 普 勒 定位 方程 ， 


a b c 


可 以 得 到 惯性 坐标 系 下 的 坐标 值 (x, y, z)。 通 过 已 有 
的 DEM 图 像 (2c) 对 其 进行 像素 偏 移 量 修正 ， 最 后 得 
到 地 理 坐 标 下 的 SAR 图 像 (2d)。 

3) 后 向 散射 系数 提取 。 在 地 面 模型 中 ， 当 
c^(0) <2 时 , 后 向 散射 系数 模型 可 以 表达 为 : 

o° (0) = K cos (1) 

式 中 : 0 为 入 射 角 ，c?(0) 为 后 向 散射 系数 。 通 过 线 
性 转换 后 向 散射 强度 到 dB 计量 后 得 到 研究 区 后 向 
散射 图 像 (2e)。 


e 


2 Sentinel-1A 图 像 提 取 后 向 散射 系数 过 程 
Fig.2 Sentinel-1A data extraction processing of backscattering coefficients 
a: 单 视 辐射 定 标 图 ; b: 多 视图 像 拼 接 ; c: SAR 坐标 下 DEM; d: SAR 地 理 编码 ; e: 后 向 散射 图 像 。a: SLC radiometric calibration; b: 
MLI mosaic; c: DEM heights in SAR geometry coordinate system; d: SAR geocoding image; e: backscatter image. 


1.2.2 ”土壤 含水 量 样本 数据 
由 于 缺乏 野外 实测 数据 ， 土壤 0-10 cm 水 分 
样本 点 数据 通过 黑龙 江 省 6 月 温差 植被 热 惯量 经 
验 模型 计算 获得 ， 其 精度 为 97.15%。 该 多 元 线性 
模型 为 : 
w=101.78—2.24AT +0.32x NDVI (2) 
式 中 : wo 为 土壤 含水 量 (%); AT 为 层 夜 温差 ('C )， 
通过 2015 ££ 6 H 21 日 同 天 的 MOD11A1 地 表 温 度 
产品 计算 获得 ; NDVI 为 归 一 化 植被 指数 (normalized 
difference vegetation index, NDVI), EH 2015 年 6 H 
16 H Landsat 8 图 像 的 近 红 外 波段 和 红 波 段 计 算 获 
得 。 根 据 地 表 覆 盖 图 ， 将 样本 点 均匀 分 布 在 不 同 土 


地 覆盖 类 型 上 (灌木 林 , 草地 及 耕地 ) 和 不 同 的 坡度 上 ， 


样本 点 数目 共计 71 个 (图 1b)。 
1.3 ”研究 方法 
1.3.1 土壤 后 向 散射 系数 提取 方法 
土壤 含水 量 、 土 壤 粗糙 度 、 植 被 的 组 合体 同时 


影响 后 向 散射 系数 ，SAR 后 向 散射 系数 与 土壤 合 水 
量 之 间 存 在 着 一 种 非 线 性 关系 "9。 植 被 种 类 与 冠 层 
结构 同样 也 与 后 向 散射 系数 有 关 忆 2)， 后 向 散射 系数 
很 难 反 演 高 植被 覆盖 度 下 的 土壤 湿度 情况 。 为 了 建 
立 土壤 后 向 散射 系数 与 土壤 合 水 量 之 间 的 关系 , 通 
过 多 模式 (多 波段 、 多 极 化 、 多 入 射 角 ) 地 表 微 波 数据 
结合 其 他 辅助 数据 (光学 遥感 、 地 表 参 数 )， 提 取 相 关 
系数 建立 模型 。 

1) 单 极 化 图 像 (VV/VH) 与 NDVI 结合 。 一 种 方式 
是 采用 合适 的 微波 植被 散射 理论 模型 ， 定量 估算 地 
表 散 射 在 总 散射 中 的 比例 。 水 云 模型 以 辐射 传输 理 
论 作 为 基础 ,考虑 地 表 植 被 覆盖 分 布 情况 , 利用 植 
被 体 散射 及 地 表 豪 减 后 散射 之 间 的 关系 结合 植被 归 
一 化 指数 NDVI 计算 出 植被 覆盖 下 的 土壤 背景 后 向 
散射 系数 。 表 达 式 所 如 下 : 


oos (0) = 07, (8) * Y^ (9)ooi(O) (3) 
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Tg (0) = Acos(8)[1— y^ (8)] (4) 


Y^ (0) = exp[-2.Bm, / cos(0)] (5) 
式 中 : ci(O) 为 土壤 背景 后 向 散射 系数 ，cves(O) 是 


植被 散射 系数 ; yy”(9) 指 植被 双 层 衰减 因子 ( 透 过 率 )， 


OKANA; A. B 的 值 取 决 于 植被 类 型 及 电磁 波 频 
率 的 经 验 常数 ; m, 是 植被 含水 率 。 
研究 区 的 地 表 履 盖 类 型 为 耕地 、 草 地 及 稀 玻 灌 
木 等 ， 其 农作物 以 小 麦 为 主 。 由 于 受 地 面 耕作 条 件 
的 限制 ， 对 地 表 植 被 参数 进行 测量 难以 实施 。 这 里 
采用 Bindlish 等 ("在 不 同 地 表 覆 盖 方 式 中 水 云 模型 
参数 ( 表 2)， 选 择 综合 利用 方式 设置 4-0.0012, B= 
0.091, 对 于 植被 含水 率 (vegetation water content, VWC), 
选择 Jackson 建立 C 波段 作物 及 植被 参数 ， 由 经 验 模 
RIDERS (RU: 
VWC-1.913 ANDVI-0.321 5NDVI (6) 
式 中 : VWC 为 植被 含水 率 ; NDVI 为 归 一 化 植被 指 
数 。 依 据 上 述 公式 则 可 提取 出 VV 极 化 及 VH 极 化 
去 除 地 表 植 被 的 土壤 后 向 散射 系数 ,为 下 一 步 进行 
土壤 含水 量 反 演 提供 合理 的 参数 。 


表 2 半 经 验 模 型 中 不 同 土地 利用 方式 的 植被 参数 值 
Table 2 Vegetation parameters values of different land use 
types in the semi-empirical model 


参数 综合 方式 草地 冬小麦 牧场 
Parameter All land-uses Grassland Winter wheat Pasture 
A 0.001 2 0.000 9 0.001 8 0.001 4 
B 0.091 0.032 0.138 0.084 
a 2.12 1.87 10.6 1.29 


2) 双 极 化 图 像 (VV+VH)。 另 一 种 方式 是 基于 统 
计 方 法 利用 地 表 双 极 化 散射 特性 与 实测 SWC 建立 
关系 。 如 : 运用 线性 去 极 化 比率 PR,ww 对 土壤 含水 
量 进 行 反 演 请， 也 可 以 是 双 极 化 系数 作 减 去 极 化 
Ow Ow RE Ow, 一 oo 进行 反 演 。 本 文 利用 Sentinel- 
1A 影像 的 双 极 化 影像 ， 分 别 通过 VV. VH 极 化 进入 
通道 进行 提取 后 向 散射 值 ， 同时 计算 图 像 的 双 极 化 
后 向 散射 系数 去 极 化 比率 。 双 极 化 图 像 分 别提 取 后 
向 散射 系数 反 演 地 表 参 数 可 以 减少 异 质 性 的 影响 ， 
每 一 种 极 化 方式 所 对 应 的 地 表 特 征 可 以 用 于 区 分 地 
物 类 型 、 土 壤 条 件 、 植 被 状况 等 ， 从 而 提高 反 演 的 
精度 。 
1.3.2 ”土壤 含水 量 反 演 模型 

水 云 模 型 能 去 除 植被 冠 层 散射 、 树 干 散射 后 的 
土壤 后 向 散射 系数 ， 当 植被 超过 一 个 临界 体积 时 ， 
衰减 地 面 散 射 逐 步 变 小 。 植 被 冠 层 的 光学 厚度 同样 


影响 土壤 的 水 分 含量 ， 模 型 加 入 NDVI 作为 辅助 变 
量 对 植被 复杂 地 区 的 土壤 水 分 值 进 行 修正 。 近 年 来 ， 
在 土壤 含水 量 估 测 中 支持 向 量 回归 (support vector 
regression，SVR) 的 应 用 越 来 越 频繁 ， 它 具有 精确 估 
计 、 易 于 使 用 、 和 良好 的 内 在 泛 化 能 力 、 处 理 非 线性 
问题 能 力 等 优点 二 。 支 持 向 量 回 归 是 将 土壤 后 向 散 
射 系 数 ( oo 或 者 c4 )、 植 被 归 一 化 指数 NDVI 与 土 
壤 水 含量 之 间 的 非 线 性 问题 , 转换 为 三 者 高 纬 特征 
空间 的 线性 问题 ， 并 求 得 最 优 分 类 面 ， 其 核心 表达 
RaR: 


I (4 -axK(, y)*b 1) 
i-l 


式 中 : a 为 样本 值 ; b 为 常数 项 ; K(x, y)79 ££ 9] 
函数 4， 能 够 将 数据 映射 到 高 维 空间 ， 其 具体 表 
K(x ERE (ô » 0) (8) 
Í 28? 

式 中 : 5 为 达 率 , 即 土壤 水 分 值 跌落 到 0 的 速度 参数 。 

利用 1.3.1 中 方法 提取 土壤 后 向 散射 系数 ， 结 合 
地 面 土壤 水 分 数据 样 点 和 地 形 数 据 ， 建 立 土壤 水 分 
空间 数据 库 ; 其 次 选择 56 个 样 点 作为 训练 集 ，15 个 
样 点 作为 测试 集合 进行 SVR 建 模 ; 依据 模型 对 应 的 
参数 ， 对 预 处 理 后 的 SAR、NDVI 分 别 进行 参数 提取 ， 
单 极 化 的 方式 输入 参数 为 VV、VH、VV+NDVI 和 
VH+NDVI， 对 于 双 极 化 的 方式 输入 参数 为 PR; 
根据 支持 向 量 回归 模型 计算 出 该 区 域 的 土壤 含水 量 ， 
并 作 精 度 验证 ; 选择 最 优 的 极 化 组 合 方式 作为 研究 
区 土壤 水 含量 反 演 结果 ,最 后 结合 其 他 辅助 数据 
DEM, 区 域 的 土壤 图 和 土地 利用 图 对 结果 进行 分 析 。 
2 结果 与 分 析 
2.1 不同 极 化 方式 下 土壤 含水 量 反 演 结果 

通过 半 经 验 的 水 云 模 型 提取 土壤 后 向 散射 系数 ， 
将 支持 向 量 回归 算法 运用 到 土壤 含水 量 模型 中 ， 反 
演出 不 同 极 化 方式 下 的 土壤 水 含量 结果 图 。 并 对 结 
果 重 新 分 为 5 类 ， 分 别 为 : 0~10%、11%~15%、16%~ 
25%、26%~30%、30%~50%， 形 成 研究 区 土壤 水 合 
量 空间 分 布 图 (图 3)。 

从 图 3 可 知 , 不 同 极 化 方式 下 研究 区 土壤 合 水 
量 结果 图 差异 性 较 大 , VH 极 化 反 演 结果 (图 3b) 数 值 
整体 偏 低 , 而 VH 极 化 与 NDVI 组 合 ( 图 3d) 反 演 却 整 
体 偏 高 。VV 极 化 、VV 与 NDVI 组 合 、 双 极 化 反 演 
结果 在 空间 分 布 上 表现 较 一 致 。 从 图 3c 可 知 ， 研 究 
区 土壤 合 水 量 大 部 分 介 于 10%~25%; SWC 为 0~10% 
区 域 一 部 分 出 现在 受 人 类 活动 干扰 较 大 的 高 速 路 周 
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3 不 同 极 化 方式 下 研究 区 土壤 水 含量 (SWC) 反 演 结果 图 
Fig.3 Soil water content (SWC) inversion results of different 
polarization ways in the study area 


的 阳 面 及 坡度 较 大 区 域 ; SWC 为 11%~15% 区 域 主要 
分 布 在 西南 小 部 分 地 区 和 北部 小 部 分 地 区 的 草地 、 
闲置 的 农用 地 、 无 植被 的 田 卉 等 区 域 ， SWC 为 
16%~25% 区 域 多 为 稀 玖 的 灌木 林 ， 土 壤 含 水 量 主要 
受 植被 蔡 腾 作用 及 地 势 条 件 的 影响 。 耕地 的 SWC 多 
集中 在 25%~30%, 主要 集中 在 河 漫 滩 附 近 的 一 级 阶 
地 上 ,用 于 耕种 农作物 ,植被 覆盖 度 较 低 ,同时 受 
灌溉 条 件 的 影响 故 其 值 较 高 。 由 于 处 于 耕作 季节 ， 
土壤 水 分 含量 为 30%~50% 区 域 位 于 河 漫 滩 低 洼 的 
区 域内 ,并 存在 沼泽 地 。 总 的 来 说 ， 该 区 域 土壤 含水 
量 状况 良好 ， 有 利于 农作物 的 生长 。 
22 ”不 同 极 化 方式 下 土壤 含水 量 反 演 精度 

在 支持 向 量 回 归 算 法 下 ， 对 应 不 同 极 化 方式 下 
反 演 结果 的 拟 合 优 度 (R 习 、 均 方 根 误差 (RMSE) 列 于 
X 3, 不 同 极 化 方式 下 的 反 演 精度 验证 图 如 图 4 所 
示 。 在 单 极 化 拟 合 方式 上 , VV 单 极 化 R2 为 0.754, 
而 VH 单 极 化 R2, 73 0.521, VV 在 单一 极 化 通道 中 对 


土壤 含水 量 反 演 表现 更 加 敏感 ， 双 极 化 RA 为 
0.595, EE VH 极 化 精度 提高 0.074; 加 入 NDVI 植被 
参数 后 ， 反 演 精 度 得 到 提高 ， 其 中 最 好 的 极 化 组 合 


表 3 不 同 极 化 方式 下 土壤 含水 量 反 演 精度 表 
Table 3 Inversion precision of soil water content in different 
polarization ways 


gml 
PE o x il EMBE 
VV 0.754 4.715 
VV«*NDVI 0.842 2.071 
VH 0.521 8.252 
VH+NDVI 0.536 7.549 
VV/VH 0.595 5.695 


VV+NDVI 


* 
e 
A 
x 
x 


反 演 值 Inversion value (96) 


5 
10 15 20 25 30 35 40 45 
SWC 样 本 值 Sample value of SWC (96) 


4 不 同 极 化 方式 下 土壤 含水 量 (SWC) 验 证 精度 
Fig.4 Verification accuracy of soil water content (SWC) 
inversion in different polarization ways 
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式 为 VV A NDVI A, Rangi 达 0.842， RMSE 为 
2.071, 比 VV 极 化 精度 提高 了 0.088, RMSE 降低 
2.704。 综 上 所 述 ， 选 择 恰 当 的 后 向 散射 参数 ， 同 时 
采用 有 效 的 辅助 参数 解决 单一 极 化 的 吸收 问题 很 重 
要 。 其 验证 精度 如 图 5 所 示 ,， 拟 合 优 度 为 : Ronai’? 
Rw>Rrwm>Rwindai>Rh， 所 以 选用 精度 较 好 的 VV+ 
NDVI 组 合 反 演 结 果 作为 最 后 地 表土 壤 合 水 量 图 。 
2.3 不同 极 化 方式 的 适用 性 分 析 
2.3.1 单 极 化 方式 反 演 结果 适用 性 分 析 

3a FB, SWC 值 较 高 地 区 分 布 在 河 漫 滩 周围 ， 
其 次 是 草地 、 灌 木林 。 在 灌木 林 区 域 SWC 值 与 实测 


值 (图 6) 相 比 低 10.45%, 由 于 植被 具有 双向 反射 与 
透射 特性 ， 是 一 个 非 朗 伯 体 结构 ,土壤 后 向 散射 为 
衰减 后 地 面 散射 ， 在 高 植被 覆盖 区 不 适合 ,其 SWC 

反 演 值 明 显 偏 低 。 图 3b 中 , VH 极 化 反 演 结果 在 灌木 
林 和 耕地 地 区 反 演 值 偏 低 ， 同 时 ， 在 河 漫 滩 耕 地 区 
域 的 土壤 含水 量 明 显 偏 低 ， 比 实际 值 平均 低 9.37%。 
根据 土壤 后 向 散射 值 与 样本 点 土壤 含水 量 之 间 的 线 
性 关系 ， 生 成 不 同 极 化 方式 下 土壤 合 水 量 的 敏感 性 
图 (图 5)。 因 为 交叉 极 化 VH 的 回 波 强 度 比 同 极 化 
VV 15, 故 VH 敏感 性 比 VV 极 化 小 , 在 一 般 情 况 下 
应 避免 选择 VH 极 化 图 像 反 演 土壤 含水 量 。 


土壤 含水 量 Soil water content (96) 
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5 VV(a)、VH(b) 极 化 对 土壤 含水 量 的 敏感 性 


Fig. 5 Sensitivity of radar signal in both VV (a) and VH (b) polarization to soil water content 


基于 VV 极 化 不 适合 在 高 植被 覆盖 地 区 使 用 ， 
VH 极 化 在 耕地 反 演 值 偏 低 ， 采 用 单 极 化 图 像 不 能 
综合 反映 地 表 履 盖 复 杂 地 区 的 土壤 含水 量 状况 。 在 
单 极 化 中 加 入 辅助 变量 NDVI 后 , 图 3d ARED, 
VH 与 NDVI 结 合 反 演 值 整体 偏 高 ,尤其 是 在 草地 区 
域 明 显 比 其 他 区 域 值 更 高 ， 比 实际 值 高 6.62%。 虽然 
NDVI 的 加 入 提高 了 低 植 被 区 域 对 水 分 的 敏感 性 ， 
但 是 在 高 植被 区 的 水 分 敏感 性 并 不 高 ; 而 VV 与 


NDVI 结合 反 演 结果 在 不 同 地 物 条 件 下 更 加 接近 真 
35r 
m 草地 Grassland 
30 


N 
Un 


N 
© 


土壤 水 分 含量 Soil water content (96) 


0 
样本 值 Sample value VV 


o $i Arable land 


VV/VH 


实 值 ， 对 于 地 表 辐 射 贡献 受到 表面 作用 支配 , NDVI 
修正 了 VV 垂直 极 化 在 植被 覆盖 区 对 水 分 的 灵敏度， 
这 一 组 合 整体 反 演 过 程 表现 良好 。 

土壤 合 水 量 与 地 表 覆 盖 方 式 有 关 ， 其 合 水 量 从 
高 到 低 为 SWCyyw>SWCy4.>SWCwj 。 适 用 于 草地 区 
域 的 极 化 方式 是 VV 极 化 .VV 与 NDVI 的 组 合 模式 ， 
其 与 真实 值 分 别 相差 2.05%、0.95%( 图 6); 耕地 区 域 
多 集中 在 河 漫 滩 附近 ， 受 灌溉 条 件 影响 ， 土 壤 合 水 
量 较 高 ， 属 于 低 植 被 覆盖 地 区 ,VV 极 化 、VYV 5 


日 灌木 Shurd 


VH+NDVI 


VH VV+NDVI 


6 不 同 土地 利用 方式 下 土壤 含水 量 反 演 结果 的 平均 值 


Fig.6 Averages of inversion results of soil water content of different land use types 
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NDVI 的 组 合 模式 同样 适用 于 该 地 区 ; 灌木 区 属于 
中 等 植被 覆盖 度 区 域 ， 在 加 入 NDVI 参数 后 单 极 化 
的 反 演 方式 更 能 够 将 植被 特性 体现 出 来 ， 而 对 比 图 
6 发 现 VV 5& NDVI £B & Eb VH 5 NDVI £8 & 88 f fc 
势 。 综 上 所 述 ， 通 过 单 极 化 结果 对 比 ， 反 演 适 应 性 对 
比 发 现 , VV 与 NDVI 组 合 更 能 综合 反 演 出 该 区 域 的 
土壤 含水 量 。 

2.3.2 ” 双 极 化 方式 反 演 结果 适用 性 分 析 

由 于 VV 和 VH 极 化 方式 彼此 强烈 的 相关 ， 在 算 
法 中 提供 相同 的 信息 ,所 以 选择 双 极 化 后 向 散射 系 
数 比率 ( PR ) 作 为 反 演 参数 能 够 尽 可 能 地 解决 数 
据 之 间 的 相关 性 。 3e 结果 显示 ,去 极 化 比率 
VV/VH 反 演 结果 精度 在 低洼 及 草地 区 域 比 VV. VH 
反 演 结果 精度 低 ,与 真实 值 相差 3.93%( 图 60)。 在 灌 
木 及 耕地 地 区 同 单 极 化 相似 表现 并 不 理想 ， 差 值 分 
别 为 11.01%、9.61%( 图 6)。 经 做 商 处 理 的 去 极 化 比 
率 反 而 造成 反 演 结果 精度 下 降 。VV/VH 组 合 降低 了 
反 演 结果 数值 ， 这 是 数据 的 相关 性 造成 的 ， 说 明 做 
商 的 方法 并 不 能 完全 解决 双 极 化 数据 之 间 的 相关 
性 。VV/VH 双 极 化 反 演 对 于 植被 覆盖 的 区 域 并 不 适 
用 ， 而 应 用 于 低 植被 覆盖 区 或 者 是 裸 土 区 域 等 的 土 
壤 合 水 量 估 测 相 对 较 合 理 。 

综 上 所 述 , 在 不 同 极 化 方式 下 进行 土壤 含水 量 
反 演 ， 对 上 比 发 现 VV/VH 双 极 化 方式 , VV 极 化 能 够 适 
用 于 裸露 和 低 植被 覆盖 地 区 的 土壤 含水 量 反 演 ， 而 
对 于 中 等 植被 覆盖 度 地 区 vv 5 NDVI 组 合 比 单一 
的 极 化 更 加 适用 。 针 对 整个 区 域 环境 的 复杂 性 ， 选 
用 能 够 综合 反映 土壤 含水 量 的 模式 是 VV 与 NDVI 


组 合 。 


3 ”结论 与 讨论 

基于 欧 空 局 Sentinel-1A 卫星 影像 ， 运 用 支持 向 
量 回归 算法 ， 对 不 同 极 化 方式 下 的 土壤 含水 量 反 演 
方法 进行 对 比分 析 ， 并 对 北 安 一 黑河 高 速 公 路 治 线 
地 区 进行 了 土壤 合 水 量 反 演 , 得 出 如 下 结论 : VV 极 
化 图 像 对 裸露 地 区 反 演 更 有 具有 优势 ,但 该 极 化 信号 
在 植被 层 会 被 强烈 吸收 ， 在 中 等 植被 覆盖 度 地 区 并 
不 适用 ; VH 极 化 反 演 效果 不 及 VV 极 化 有 效 ; 双 极 
化 去 极 化 比率 在 低 植被 覆盖 地 区 更 能 发 挥 优势 ， 如 
草地 、 牧 场 等 区 域 ,但 不 能 应 用 到 高 植被 覆盖 地 区 :; 
而 VV 与 NDVI 组 合 更 能 够 综合 反 演出 该 区 域 复杂 
地 表 的 土壤 合 水 量 , RIBA 1 号 C 波段 SAR 影像 
的 极 化 特性 优势 。 

不 同 极 化 组 合 方式 ,在 微波 遥感 反 演 植被 覆盖 
区 土壤 水 分 时 的 影响 不 同 。VYV 极 化 方式 能 将 垂直 信 


息 反 映 良 好 ,在 植被 覆盖 区 域 上 提取 的 土壤 后 向 散 
射 系数 一 般 值 偏 低 ， 反 演 值 也 会 偏 低 。 引 入 NDVI 
作为 修正 植被 影响 的 部 分 ,可 将 植被 -土壤 之 间 的 非 
线性 关系 表达 得 更 为 详细 、 精 确 。VH 交叉 极 化 包含 
垂直 及 水 平方 向 的 信息 ,水 平方 向 的 信息 对 于 裸露 
地 表 更 适用 ,这 与 Bindlish 等 (小 、Pasolli 等 253 针 对 
不 同 极 化 SAR 反 演 土壤 含水 量 研究 结果 一 致 。 支 持 
向 量 回归 算法 能 够 精确 反 演 土壤 含水 量 参数 之 间 的 
非 线 性 关系 , 但 是 算法 对 于 大 规模 的 反 演 来 说 需要 
较 大 的 运算 内 存 和 时 间 P "在 光学 植被 指数 的 选 
取 中 只 用 了 光学 厚度 这 一 方面 ， 并 未 考虑 其 他 植被 
旨 数 及 其 他 相关 因子 。 因 此 在 未 来 的 研究 中 , 应 该 
结合 区 域 特 点 综合 考虑 不 同 植被 指数 、SAR 不 同 入 
射 角 及 气象 因子 对 地 表 水 分 的 影响 ,进一步 完善 地 
表土 壤 合 水 量 模型 。 
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